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We derive first next-to-next-to-leading logarithmic resummations for jet-veto efficiencies in Higgs 
and Z-boson production at hadron colliders. Matching with next-to-next-to-leading order results 
allows us to provide a range of phenomenological predictions for the LHC, including cross-section 
results, detailed uncertainty estimates and comparisons to current widely-used tools. 



In searches for new physics at hadron colliders such as 
the Tevatron and CERN's Large Hadron Collider (LHC), 
in order to select signal events and reduce backgrounds, 
events are often classified according to the number of 
hadronic jets — collimated bunches of energetic hadrons 
in the final state. A classic example is the search 
for Higgs production via gluon fusion with a subsequent 
decay to W+W~ [J 0. A severe background comes 
from ft production, whose decay products also include 
a W + W~ pair. However, this background can be sep- 
arated from the signal because its W + W~ pair usually 
comes together with hard jets, since in each top decay 
the W is accompanied by an energetic (6) quark. 

Relative to classifications based on objects such as lep- 
tons (used e.g. to identify the W decays), one of the dif- 
ficulties of hadronic jets is that they may originate not 
just from the decay of a heavy particle, but also as Quan- 
tum Chromodynamic (QCD) radiation. This is the case 
in our example, where the incoming gluons that fuse to 
produce the Higgs quite often radiate additional partons. 
Consequently, while vetoing the presence of jets elimi- 
nates much of the ft background, it also removes some 
fraction of signal events. To fully interpret the search 
results, including measuring Higgs couplings, it is cru- 
cial to be able to predict the fraction of the signal that 
survives the jet veto, which depends for example on the 
transverse momentum threshold pt.veto used to identify 
vetoed jets. 

One way to evaluate jet-veto efficiencies is to use a 
fixed-order perturbative expansion in the strong cou- 
pling a s , notably to next- next-to- leading order (NNLO), 
as in the Higgs-boson production calculations of Q- 
5]. Such calculations however become unreliable for 
Pt.veto *C M, with M the boson mass, since large terms 
a™L 2n appear [L = ln(M/p t ,veto)) in the cross-section 
to all orders in the coupling constant. These enhanced 
classes of terms can, however, be resummed to all or- 
ders in the coupling, often involving a functional form 
exp(Lgx(a s L) + g 2 (a s L) + a s g 3 (a s L)/Tr + . . .). 

There exist next-to-next-to-leading logarithmic 
(NNLL) resummations, involving the g n (a s L) functions 
up to and including g 3 , for a number of quantities 
that are more inclusive than a jet veto: e.g. a Higgs 



or vector-boson transverse momentum (m9(, the beam 
thrust [l(J, and related observables [lllTlil]. To obtain 
estimates for jet vetoes, some of these calculations 
have been compared to or used to reweight [T3, Il3l4lt| 
parton-shower predictions [lT], Gil matched to NLO 
results pjl, Hoj]. However, with reweighting, neither the 
NNLO nor NNLL accuracy of the original calculation 
carry through to the jet veto prediction. 

Recently there has been progress towards NNLL cal- 
culations of the jet veto efficiency itself. Full NLL results 
and some NNLL ingredients for Higgs and vector-boson 
were provided in [2l| . Ref. [22| used these and other in- 
gredients in the soft-collinear effective theory framework 
to consider resummation for the Higgs-boson case beyond 
NLL accuracy. In this letter we show how to use the re- 
sults of [2l| together with those from boson p t resumma- 
tions [6-9] to obtain full NNLL accuracy. We also exam- 
ine the phenomenological impact of our results, including 
a matching to NNLO predictions. Given the ubiquity of 
jet cuts in hadron-collider analyses, the understanding 
gained from our analysis has a potentially wide range of 
applications. 

The core of boson transverse-momentum (pf ) resum- 
mations lies in the fact that soft, collinear emissions at 
disparate rapidities are effectively emitted independently. 
Summing over all independent emissions, one obtains the 
differential boson p t cross section 
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where <tq is the leading-order total cross section, 
\dki\M 2 (ki) is the phase-space and matrix-element for 
emitting a soft, collinear gluon of momentum ki, while 
the exponential factors and b integral encode in a fac- 
torised form the constraint relating the boson p t and 
those of individual emissions 5 2 {pf — Y17=i ^a) [23]. The 
— 1 term in the round brackets arises because, by uni- 
tarity, virtual corrections come with a weight opposite 
to that of real emissions, but don't contribute to the pf 
sum. 

To relate Eq. jjl to a cross section with a jet-veto, 
let us first make two simplifying assumptions: that the 
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independent-emission picture is exact and that a jet al- 
gorithm clusters each emission into a separate jet. The 
resummation for the cross section for the highest jet pt 
to be below some threshold pf, considering jets at all ra- 
pidities, is then equivalent to requiring all emissions to 
be below that threshold: 
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with the same universal matrix element M 2 (ki) entering 
Eqs. © and ©. 

Eq. ((2|) is clearly an oversimplification. Firstly, even 
within the independent emission picture, two emissions 
close in rapidity y and azimuth (j) can be clustered to- 
gether into a single jet. Let us introduce a function 
J(k\,k 2 ) that is 1 if fci and k 2 are clustered together 
and otherwise. Concentrating on the 2-emission con- 
tribution to Eq. (JSJ), one sees that clustering leads to a 
correction given by the difference between the veto with 
and without clustering: 

j*h»t aQ = ^J [dk 1 ][dk 2 ]M 2 (k l )M 2 {k2) x J(h,k 2 ) 

(9(p t J ~ k tA2 ) - Q(p- t ' - k ht )Q(pf - k 2 ,t)) ■ (3) 

where k\ 2 — k\ + k 2 (throughout, we assume standard 
_E-scheme recombination, which adds 4- vectors). This 
contribution has a logarithmic structure a 2 L, i.e. NNLL, 
with each emission leading to a power of a s , while the 
L factor comes from the integral over allowed rapidities 
{\y\ < ln(M/p t>veto )). 

For more than two emissions, two situations are possi- 
ble: (1) three or more emissions are close in rapidity, giv- 
ing extra powers of a s without extra log-enhancements 
(N 3 LL and beyond); (2) any number of extra emissions 
are far in rapidity, each giving a factor a s L, i.e. also 
NNLL. The latter contribution is simple because, inde- 
pendently of whether the two nearby emissions clustered, 
those that are far away must still have pu < pt, veto- 
Thus the full "clustering" correction to the independent- 
emission picture is a multiplicative factor (1 + J rclust ), as 
first derived in detail in the appendix of [21| using results 
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For the generalised- k t jet-algorithm family [25|-|29j, 
with a jet radius parameter R, we have J(ki,k 2 ) — 
Q(R 2 - (yi - 2/2 ) 2 - (0i - </>2 ) 2 )- At NNLL accuracy 
Eq. (^evaluates to J" clust = 4a2(p t ^ vcto )C* J L/ clust (i?)/ 7 r 2 
with [2l| 
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for R < it; C is Cf = 
incoming quarks (e.g. qq 
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, . . = 3 respectively for 
— > Z) or incoming gluons (e.g. 



Next, we address the issue that gluons are not all emit- 
ted independently. This is accounted for in Eq. (p} be- 
cause, to order a 2 , the M 2 (k) quantity that appears there 
should be understood as an effective matrix element 

[dk}M 2 (k] = [dk](M 2 (k) + Mt loop (k)) 

d 2 k t [dk a }[dk b }M 2 OIYcl (k a , k b )5 2 (k t . ab - k t ) (5) 



where M 2 (k) is the pure O (a s ) matrix element, 
M 2 orrcl (k a , k b ) the correlated part of the matrix element 
for emission of two soft-collinear gluons or a quark- 
antiquark pair, including relevant symmetry factors, and 
M 2 _ loop the corresponding part of the a 2 1-loop matrix 
element. The separation into correlated and independent 
emissions is well defined because of the different colour 
factors that accompany them in the generic case (30l - [33l | . 
The (5-function in Eq. ([5]) extracts two-parton configura- 
tions with the same total pt as the f-gluon configurations. 

For a jet veto, part of the result from the effective ma- 
trix element carries through: when two correlated emis- 
sions are clustered into a single jet, it is their sum, k t , a b, 
that determines the jet transverse momentum. There- 
fore the same effective matrix element can be used in 
Eq. @, as long as one includes an additional correction 
to account for configurations where the two emissions are 
clustered in separate jets: 

J-° rre V = a f[dk a ][dk b ]M 2 orIel {k a ,k b )x(l-J(k a ,k b )) 
(Q(p J t - k ta )Q(p J t ~ hb) ~ Q{p J t - h,ub)) ■ (6) 
At NNLL, F coml = 4a 2 (p t ,vcto)C J L/ corrol (i?)/7r 2 with 

™i,m ,'(-131 + 127r 2 + 1321n2)CU 
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+ (23 = t ln2)n ^ 4+0.61^-0.015^+0 (R 2 ) , 
72 J R 

(7) 

for generalised-fc/ algorithms, in the limit of small R. 
Ref. 21] includes a numerical result for all R < 3.5 and 
analytical terms up to i? 6 , used in the rest of this article. 
It did not, however, make the relation with the boson p t 
resummation. 

All remaining contributions to a NNLL resummation, 
such as the 3-loop cusp anomalous dimension or a mul- 
tiplicative C\a s term are either purely virtual, so inde- 
pendent of the precise observable, or involve at most a 
single real emission, so can be taken from the boson p t re- 
summations Q Thus the full NNLL resummed cross 



1 For generic processes, subtleties can arise with spin-correlation 
effects 34j. These are simpler for jet vetoes, which don't correlate 
distinct collincar regions. 
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section for the jet-veto is given by: 
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(8) 

where the coefficient functions "H^ 1 -* and Ci^, and re- 
summation functions gun? and (73 are as derived for the 
boson pf resummation [6|,LZ|,|9[ (reproduced for complete- 
ness in the supplemental material to this letter [35|, to- 
gether with further discussion on the connection to boson 
pf resummation). The results are expressed in terms of 
L = ln(Q/p t!VC to), ol s = a s (nFi)] the resummation, renor- 
malisation and factorisation scales Q,fiR and fi F are to 
be chosen of order of M. 

A form similar to Eq. ([8]) was derived independently in 
22] for Higgs production, also using ingredients from |2l| . 
It differs however at NNLL in that the combination 
£ jciust _|_ jcorroi j g acconi panied by an extra — (,?,Ca- 
Ref. [22| had used a NNLL analysis of the R — > 00 limit 
to relate jet and boson-pt resummations. A subtlety of 
this limit is that one must then account for a N 3 LL a 2 R 
term, which for R > In M/p t is promoted to an additional 
NNLL a 2 \nM/p t contribution f35| . 

One check of Eq. (J5J is to expand it 
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in powers of 

(Pt) = a?Er=o S NNLL,nbt)> and compare 
(pt)/d\np t to the NLO Higgs+l jet predic- 

from MCFM dT, { 2 J) (p t )/d\np t . NNLL 
resummation implies control of terms a 2 L 3 . . . a 2 (con- 
stant terms) in this quantity and so the difference be- 
tween MCFM and the 2nd order expansion of the resum- 
mation should vanish for large L. This is what we find 
within reasonable precision. The precision of the test 
can be increased if one considers the O (a 2 ) difference 
between the jet and boson-p t resummations, which has 
fewer logarithms and so is numerically easier to deter- 
mine in MCFM. It is predicted to be 
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(f clust (R) + f com \R) + (3C) . (9) 



This is compared to MCFM's LO H+2-jet result in the 
upper panel of Fig. [TJ There is excellent agreement at 
small pt, for each of three R values. The result of [22[ 
(BN, only for R — 0.5) is also shown for comparison. 
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FIG. 1. Upper panel: 2nd order difference between jet and 
Higgs-boson lnpt differential distributions, showing the coef- 
ficient of 4ttsCA(To/T 2 as determined with MCFM and pre- 
dicted in Eq. (|9l). for three R values. We also show the pre- 
diction from [23 (BN). Lower panel: differences at O (o^oo) 
between jet and boson lnpt differential distributions, with the 
expected a 3 8 aoL 2 term subtracted (denoted by a subscript ij n ), 
showing the MCFM H+2-jet NLO result compared to our 
NNLL prediction for the o? s <tqL term. 



The above test can be extended one order further 
by examining the order o^oo difference between the jet 
and boson p t differential distributions. The comparison 
between our predictions and MCFM H+2-jet NLO re- 
sults [13, Ell is given in Fig. [T] (lower panel) , for each of 
three R values. To facilitate visual interpretation of the 
results, the expected a^.aoL 2 term has been subtracted. 
The residual a^a^L term is clearly visible in the MCFM 
results and, within the fluctuations, coincides well with 
our predictions, providing a good degree of corroborating 
evidence for the correctness of our results beyond order 
a 2 s cr . 

To illustrate the phenomenological implications of 
our work, we examine the jet veto efficiency e(pt) = 
^ J \Pt) I '° 'tot, where atot is the total cross section, known 
up to NNLO [44U491 ]. We combine ("match") the resum- 
mation with fixed-order predictions, available from fully 
differential NNLO boson-production calculations 0, H, 
l50l . l5lT ] or NLO boson+jet calculations [H, [52j imple- 
mented in MCFM [53|. We use three matching schemes, 
denoted a, b and c, straightforward extensions (3f| of 
those used at NLL in 1211] . 



4 




FIG. 2. Comparison of NNLO, NLL+NNLO and NNLL+NNLO results for jet-veto efficiencies for Higgs (left) and Z-boson 
(right) production at the 8 TeV LHC. The Higgs plot includes the result from a POWHEG (revision 1683) [2(1 [H} plus Pythia 
(6.426) fl7l . El simulation in which the Higgs-boson pt distribution was reweighted to match the NNLL+NNLO prediction 
from HqT 2i)Q as in [2l[. The lower panels show results normalised to the central NNLL+NNLO efficiencies. 



Our central predictions have [Ir = /if = Q = M/2 and 
scheme a matching, with MSTW2008NNLO PDFs [HI]. 
We use the anti-fc t [1^| jet-algorithm with R = 0.5, as 
implemented in Fast Jet [55|. For the Higgs case we use 
the large mt op approximation and ignore bb fusion and 
&'s in the gg — > H loops (corrections beyond this approx- 
imation have a relevant impact fl6l [56j). To determine 
uncertainties we vary /ir and /if by a factor of two in 
either direction, requiring 1/2 < [j,r/{j,f < 2. Maintain- 
ing central (J.r,f values, we also vary Q by a factor of 
two and change to matching schemes b and c. Our final 
uncertainty band is the envelope of these variations. In 
the fixed-order results, the band is just the envelope of 
fJ-R.F variations. 

The results for the jet- veto efficiency in Higgs and Z- 
boson production are shown in Fig. [2] for 8 TeV LHC 
collisions. Compared to pure NNLO results, the cen- 
tral value is slightly higher and for Higgs production, the 
uncertainties reduced, especially for lower pt.veto values. 
Compared to NNLO+NLL results [2l|, the central values 
are higher, sometimes close to edge of the NNLO+NLL 
bands; since the NNLO+NLL results used the same ap- 
proach for estimating the uncertainties, this suggests that 
the approach is not unduly conservative. In the Higgs 
case, the NNLO+NNLL uncertainty band is not particu- 
larly smaller than the NNLO+NLL one. This should not 
be a surprise, since [2l[ highlighted the existence of pos- 
sible substantial corrections beyond NNLL and beyond 
NNLO. For the Higgs c ase, we also show a prediction 
from POWHEG HJgjJ interfaced to Pythia 6.4 [l]} at 
parton level (Perugia 2011 shower tune [4J|), reweighted 



to describe the NNLL+NNLO Higgs-boson p t distribu- 
tion from HqT (v2.0) Q, as used by the LHC experi- 
ments. Though reweighting fails to provide NNLO or 
NNLL accuracy for the jet veto, for pt,veto scales of prac- 
tical relevance, the result agrees well with our central 
prediction. It is however harder to reliably estimate un- 
certainties in reweighting approaches than in direct cal- 
culations. 

Finally, we provide central results and uncertainties 
for the jet-veto efficiencies and 0-jet cross sections (in 
pb) with cuts (in GeV) like those used by ATLAS and 
CMS, and also for a larger R value: 
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Interestingly, the R = 1 results have reduced upper un- 
certainties, due perhaps to the smaller value of the NNLL 
f(R) correction (a large f(R) introduces significant Q- 
scale dependence). The above results are without a ra- 
pidity cut on the jets; the rapidity cuts used by ATLAS 
and CMS lead only to small, < 1%, differences pll ]. 

For the 0-jet cross sections above, we used total 
cross sections at 7 TeV and 8 TeV of 15.3!1'2> pb and 
19.511 i pb respectively [57], [H[ (based on results in- 
cluding 45-49]) and took their scale uncertainties to be 
uncorrelated with those of the efficiencies. Symmetris- 
ing uncertainties, we find correlation coefficients between 
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the O-jet and > 1-jet cross sections of —0.43 (—0.50) for 
R = 0.4 (R = 0.5), using the covariance matrix in |35| . 

Code to perform the resummations and matchings will 
be made available shortly. 
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Note added: as our manuscript was being finalised, 
Ref. [67jj appeared. It claims issues in NNLL resumma- 
tions of jet vetoes, however does not address the all-order 
derivation of the NNLL /^-dependent terms in [2l|. Its 
claim is further challenged by the cef numerical check in 
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SUPPLEMENTAL MATERIAL 



We here provide material that completes the discussion of the letter, including more explicit formulae, some deriva- 
tions and supplementary figures. 

1. Explicit resummation formulae 

In the present section we report the explicit expressions for the resummation functions g\ 1 gi and 53 computed 
in 0,0]) as functions of A = a s (3oL, with L = \n(Q/pt). a s denotes a s (liR) unless otherwise stated, and Q is the 
resummation scale (see main text) 

, x A^> 2A + ln(l - 2A) , , 

5i(A) = — A- '-, (10a) 

Ttpo LA 

1 , / M 2 m \ 2A+ (1 - 2A)ln(l - 2A) 

■ 92(A) = 2^ ^ - 2A > { Ai ) ln W + B ) - I - 2A 1 

fa ln(l-2A)((2A-l)ln(l-2A)-2)-4A 1 (2A(1 - ln(l - 2A)) + ln(l - 2A)) 

+ 4tt/33 l - 2A ~~ 2tt/3 1 - 2A n Q 2 

(10b) 

„ w . h - + B ,„) ( _ _A_ _ I + ^HilJ^) ) _ J__A_ ( A » ,„ £ + W 

A< 3 > A 2 (2) / gi 2A(3A-l) + (4A-l)ln(l-2A) 1 A 2 /i| N 

" 47r 2 /3 2 (1 - 2A) 2 + ^4^3 (1 - 2A) 2 ^ (1 - 2A) 2 n 

(1) / A(/?o/3 2 (l-3A)+/3 2 A) (1 - 2A) ln(l - 2A) (/3 /3 2 (l - 2A) + 2(3 2 \) (3 2 (1 - 4A) ln 2 (l - 2A) 
1, /3 4 (1-2A) 2 + 2/3 4 (l-2A) 2 + 4/3 4 (1 - 2A) 2 
in 2 # - ^ (2A(1-2A) + (1-4^)^1-2^) ^ ^ (iQc) 



(1-2A) 2 Q 2 2/3 2 (1-2A) 2 Q 2 /' 

where, for Higgs, A^> = 2C A and #W = -4tt/3o, while for Drell-Yan, A^ = 2C F and #W = -3C F . The remaining 
coefficients can be expressed in a unique way as 0, l59l |60] | : 

A( 2 ) = A«4 1 2 IW , AW^^W^w + tt/SoC^), B< 2 > = -2 7 < 2 ) + 2^ CC 2 , (11) 

= 11CU - 2ry = 17C 2 - 5C A nj - 3C F n f 

P0 12tt ' Pl 24tt 2 ' 1 ' 



2857C% + (54Cf - 615C F CU - 1415C^)fi/ + (66C F + 79C* A ) 



2 



ft- -—-- («) 

in terms of the Casimir C = Ca for Higgs and C = C F for Drell-Yan, and of the well known constants 

^^(^SCsVf^, (14a) 



K {1) 


= C A | 


f67 
K 18~ 


4)" 


5 


^CMW 


= c\[ 


'245 
v 24 ' 


-fc- 





Here 7^ 2 - 1 [6l], [62| are the coefficients of the 5(1 — z) term in the NLO splitting functions PW. For Higgs production 
we have 



7 (2) = C\ Q + 3C 3 ) - ^Cjrn/ - j^n, . 



(15) 



whilst for the Drell-Yan process 



-C r (|-i- + 6ft) + C F C. 4 g + ^ - 3< 3 - C,„, 1 + i.) . (16) 
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W W = tf« _ B W + ±_ In ^_ m ^_ + q 2^ In ^ , (17b) 



We finally report the expressions for the collinear coefhcient function C^\z) and the hard virtual term in 
eq. © 

C\f (z) = -P^{z) - 5^6(1 - z)C^ + (z) ln J, (17a) 

aw m 2 \ m 2 - 2 

T~ Q 2 ) Q 2 " ' H ^"^" i M 2 ' 

where q is the a s power of the LO cross section (q = 2 for Higgs production and q = for Drell-Yan) . The coefficient 
encodes the pure hard virtual correction to the leading order process, it is given by 

Higgs : H<U = C A \5 + ^ 2 ) - 3C F , (18a) 
Drell - Yan : H {1) = C F ^-8 + ^tt 2 ^ . (18b) 



Finally, P^ (z) is the 0(e) term of the LO splitting function P\ (z): 

P^(z) = -C F (l-z), (19a) 

P^(z) = -C F z, (19b) 

P^(z) = -z(l-z), (19c) 

P$' e (z) = 0. (19d) 

2. Full matching formulae 

We start by recalling the three prescriptions discussed in [21] for defining the jet-veto efficiency at NNLO accuracy. 
In this section, in order to simplify the notation, we will refer to the integrated jet-veto distribution T,^ as E. 

( a ), s E (pt,veto) + Ei(j3 t ,veto) + E 2 (pt,vcto) > 

; (pt,vcto) = ; ; , (20a) 

£ w {pt vcto) = E (p t , veto ) + E 1 (p t , vcto ) + E( J) (p t , vcto ) (2Qb) 

(c)f \ i I 1 (Pt.vctoJ CTi-(j) 
£ V ; (Pt,vcto) = H 2 S 1 (Pt.vcto) 



^2 (Pt,veto) 



(TO CTg (TO 

with 

,(J) 



(20c) 



Ei(p t ,veto) = <y% + E> (Pt.veto) (21) 

being the C(a*)— th correction relative to the Born cross section where 

£fWeto) = -r *t^fH (22) 

can be determined from MCFM, while <Xj is the i th order contribution to the total cross section (cf. [44l - |49l ]). The 
above three prescriptions differ by terms 0(a^), which are beyond the control of current fixed-order calculations. 

The jet-veto matched efficiency should tend to one and the differential distribution should vanish at the maximum 
allowed transverse momentum p™veto 



de 

dpt,veto 



e (p££o) = 1, - (K" a c x to) = 0. (23) 



2 Often in the literature, the hard coefficient ff W is considered as different convention for the resummation coefficient B< 2 ) which 

part of the 5(1 — z) term in the coefficient function c\- (z), so it will differ by an amount 2-wf)oH i - 1 '> from what reported here, 

comes with a factor 1/(1 — a s f3oL) in eq. J5}. This results in a 
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To fulfil such requirements, we modify the resummed logarithms as follows 



L ->L = { 1 



Pt.vcto \ 1 



Ft, veto/ P 



-In 



Q 



Pt.veto 



Q V 



„max 
-Pt,veto 



(24) 



is necessary to fulfil 



where p is some integer power. By default we choose p = 5 [2l|. The factor ^i — pAl° to 

eq. (f2"3")) but it is largely irrelevant in practice since p™vet * s mucri larger than the typical values of the jet transverse 
momentum veto (in practice, we set Pt.max — oo). We introduce three multiplicative matching schemes jr33l ]. each of 
them corresponding to one of the three efficiency definitions (|20al) . j2Qb|, (|20cp . To simplify the notation, we split 
the luminosity factor in the square brackets of Eq. (|SJ) into two terms D-°'{L) and C^(L), which start at order a° 
and a] respectively, 



id J 



x 2 ,e ii F ), 



(25) 



£(1)( ^ ) = £5Z j ' dxidx 2 5(x lX 2s - M 2 ) f i (x u e- i ti F ^f j (x2,e- L n F )uW 



+ n 7T- a F Ef /' (*)/<(-. fj(x2,e- £ ^ F ) +{(x 1} i) o (* 2 ,j)} 

l-2a s poL k \J Xl z \z / \ 

The first of the three matching schemes then reads 



(26) 



^matched (Pti veto) 



1 SNNLL(Pt,vcto) 

^1 + £( 1 )(L)/£(°)(L) 



00 



1 + + S(1)(Pt '™ to) - S NNLI>t,vcto) 



-s (2) bt,vc to )-s 



fo ^ , ( £(1) (0) ^NLL^vctoA / (1) 
'NNLL lPt,vetoJ + I £(0)(Q) — I U>t,vetoJ ~ 



^NNLL,(Pt.vcto 



and the corresponding jet- veto efficiency is 



O) / \ _ ^matched (Pt, veto) 

e matched^t,vctoj - , } 



(27) 



(28) 



(-max \ 
^matched V-Ft , veto I 

The second scheme can be derived from the previous one by replacing Y,( 2 \p t>veto ) with E^ 2 ) (pt,veto)- For the vetoed 
cross section we get 

( \_ 1 ^NNLLyPt.vcto) f / £ (1) (L) \ (1) (i) 

S ma tchcd(Pt,vcto) - ^ + [ CT0 ^ + #5)(£) J + ^ (Pt ' VCto) _ S NNLI>t,vcto) 



4-y( 2 )(n \ V (2) r« ^ ■ ( ^NNLL (PtjVeto) | Ml), x , x 

+ ^ v HPt.vctoj - ^NNLLiPt.vctoJ + I £(0)77j\ ~ I ^ U>t,vetoJ ~ ^NNLLU^vctoJ 



while for its efficiency 



e matchod(Pt. vct o) — 



^matched IPt.vcto 

•y( b ) /.max \ 

^matched l-Ft , veto I 



(29) 



(30) 



Finally, the third matching scheme is directly formulated for the efficiency resulting in 



£ (o) 
matched 



(Pt,veto) — 



f ENNLL,(pt,vcto) 

4 l + £( 1 )(Z)/£(°)(i) 



1 + £(°)(L) j + S(1) ^^ vcto ) ~~ ^NNLL (Pt,veto) 



+ E( 2 )( Pt , voto ) - ^£W(pt )W to) - 4n] 

/Tn 



"'NNLL (Pt,vetoJ 



^(O) S&UCPt.vrto) ' 

£(°)(0) cto 



(E (1) (p t ,veto) - E^ LL (p t ,veto) 



(31) 
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3. Details of relation between jet and boson-pt resummations 



This section collects a number of results to help relate jet and boson p t resummations. Firstly we demonstrate that 
the g n (a s L) from a boson p t resummation can be directly carried over to jet p t resummation for n < 3. Then we 
obtain a form for the boson p t resummation that is suitable for expansion and c omp arison with fixed-order results. 
Finally we consider the large R limit of the jet-pt resummation, which was used in [22j to attempt to obtain a relation 
between jet and boson p t resummations. 



a. Relating g n (a„L) between boson and jet resummations 



One of the main ingredients of our results are the g n (a s L) functions that are used in boson p t resummations. These 
stem from the rightmost integral in Eq. ([1]), which involves a Fourier transformation, whereas in Eq. ([2]) we need 
related integrals but with a theta-function instead of the (exp(ib.k t ) — 1) factor. 

We start from the expression for the resummed p t distribution in eq. ([1} and concentrate on the part of the matrix- 
element in the right-hand integral that is responsible for the leading logarithms. Integrating over azimuthal angles 
we obtain: 



= cto / bdbJ a (bp t ) exp[-K(b)] , 11(b) = / [dk]M 2 (k)(l - J (bk t )) . 



d^ B \p t ) 
Ptdpt 

We wish to show that we can safely perform the replacement 

(1 - J (bh)) -> 6(fc f - bo/b) , b = 2e^ E , 
up to and including NNLL accuracy. Integrating over rapidity 71(b) has the form 



M 



dk t 



M 



K(b) = / -r~F a s In — (1 - J (bk t )) 



! as ~k t 



4C^ln M 



7T k t l- 2a s /3 In f- 



(32) 



(33) 



(34) 



To evaluate separately real and virtual contributions in eq. (|34|) . we introduce a dimensional regulator and write 



U(b) = F (a s 8 e ) 



M 



dk t 
k t 



(1 - Jo(bk t )) 



e=0 



which yields 

K(b) = R hh (b /b) + 511(b), 
where, neglecting terms suppressed by powers of l/(bM), 



RLL(b /b) = 



M 



dk t 



F a. In 



M 



\ eft* - 



e(fct - b /b) . 



and 



811(b) = F (a s d e ) 



(b/boT 



"7BE 



r(i-f) 



r(i + f) 



e=0 



F(aA) [ - 



Q Y/ + o(^) 



(35) 

(36) 

(37) 

(38a) 
(38b) 



e=0 



This gives at most a term a™ ln™ _2 (Af6/6 ), i-e. a N 3 LL term. A similar argument can be applied to contributions 
to 11(b) arising from less singular regions, giving also rise to terms that are beyond NNLL. Consequently, to NNLL 
accuracy, the same g± , 172 and g% functions can be used in both the jet and boson resummation. 
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b. Evaluation of the boson-pt integrated cross section 



To facilitate comparisons between the jet and boson p t resummations at fixed order, it is convenient to have an 
expression for the boson pt resummation whose fixed-order expansion can be straightforwardly obtained. The full 
expression for the cumulative pt cross section can be found in [y, @, Q and reads 



S (s) ( Pt ) = f°° dyJi(y) \M B \ 2 e- R( -»°/V (c^(\n(Qb/b )) + £ (1) (ln(Q6/M)) 



where 



R{b /b) = HQb/b ) 9l (a s HQb/bo)) + g 2 {a s HQb/b Q )) + —g 3 (a s \n(Qb/b )) 



(39) 



(40) 



is the full NNLL radiator. As discussed above, the resummation functions gi, g 2 and 173 are those used for the jet veto 
case. To perform the inverse Fourier transform we expand R(bo/b) and the full luminosity factor around b = bo/pt 
and neglect subleading logarithmic terms getting, at NNLL accuracy, 



C^iHQM) + C^iHQ/Pt)) + d lnpt C^(\n(Q/ Pt )) \n(y/b ) 

-R' 



-R(pt 



l-ii?"ln 2 (y/6 ) j . 



(41) 



where we have performed the change of variable y = bp t , and we have made use of R' and R", the first and second 
derivatives of R with respect to ln(Q/p t ). To order a s L, R' — Aa s Cln(Q/p t )/n. Moreover, from eq. ([25]) . we see that 
the variation of £(°'(L) reads 



d lnpt C ( -°\L) = -^Yl I d Xl dx 2 5{ Xl x 2S - M 2 ) [(P^ <E> ^ (x 1 ,e- L f i F ) f J (x 2 ,e- L t i F ) + {( Xl ,i) o (x 2 ,j)}\. (42) 



It is straightforward to show that eq. (|4"Tj) evaluates to 



£ (B) (Pt) = \M B \ 2 e- R ( pt) 



£W(ln(QM)) (l-\B!'&k\ 

+ C^(HQ/p t )) - d inpt £^(\n(Q/p t ))d R , 



r(i + f) 



(43) 



In this notation, the result for the jet-veto cross section is simply \M B \ 2 e- R{pt) (£ {0) + £^)(1 + T clust + T coml ). 
It is therefore immediate to evaluate the differences between the two formulae at any given fixed order and in 
particular to derive Eq. ([9]): making use of the fact that e~ lEH L(l — -^-)/F(l + -^-) has an expansion of the form 



1 + %R' 3 + O ( R'° ) , one sees that the only terms in the difference that survive at order a 2 s L are the J rclust and J-' 



rcorrcl 



contributions and the R"d B , term of Eq. (j43[) . with the latter giving 



1 , rfl - ^ 
a Q -R"d 2 RI e -r^iii & 

2 R r(i + f 



2 r o\ 



00 -4-§( 3 C z \n^ + 0(aiL 
ft Pt 



OialL 2 \ 



(44) 



which is the source of the Cs in Eq. ©0 



3 One point to note in evaluating the difference between the jet 
and boson pt resummations at order erf L 2 is that it is necessary 
to account also for the difference between C2 terms for the two 
resummations. One of the properties of this difference of C2 
terms is that is has Q dependence that ensures that the final 



prediction for the difference of a 3 s L 2 terms is Q-independcnt. 
To produce figure [T] the difference of C2 terms was taken from 
a numerical determination based on the MCFM leading-order 
H + 2-jet calculation. 
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c. Use of the large-R limit to relate boson and jet-p t resummations 

One natural way of relating jet and boson-p t resummations is to make the observation that for an infinite jet radius, 
all partons will be clustered into a single jet ; which will have a transverse momentum that balances exactly that of 



the boson. This approach was taken in Ref. 22] and here we examine it in detail. 

First, let us consider the properties of J^ dust and J rcorrel for large R. It is straightforward to see that ^' correI vanishes 
for large R, since in Eq. ([7]) the two partons will always clustered together, giving 1 — J(ki, k 2 ) = 0. For J 70 lust , the 
NNLL component for R > it can be evaluated in closed form and is given by 

= ( ^' - £) — ■ ^ + (? " i) ^) ■ <«> 

This has the property that it vanishes as 1/R for large R. Thus it would appear that at order a 2 L the difference 
between jet and boson-p t resummations should be given by [22| 

^^NNLL,2 (Pt) ^ S NNLL,2(Pt) ferns tf \\ 2 tins 1 f A R\ 

jj -r\ = if(R) ~ /(ooj) a^ao = f{R) u s (Ta , (46) 

which differs from the result in Eq. © (here f(R) = f COIxel (R) + f clust {R)). 

To understand the origin of this difference, it is helpful to examine the structures that lead to J^ clust vanishing for 
large R. A first observation is that for large R, J(ki, k 2 ) can be written as 

J(k 1 ,k 2 ) = @(R-\Ay\ + ^ + o(^]) , Ay = yi -y 2 , Acf> = fa - fa . (47) 



2R \R 3 _ 

Neglecting the term of order 1/R will allow us to simplify our discussion and so we will instead examine a "rapidity- 
only" jet algorithm with the clustering function 

J iap (k 1 ,k 2 ) = Q(R-\Ay\). (48) 

Let us now evaluate Eq. ([3]) with J rap . We break the problem into rapidity, transverse momentum and azimuthal 
integrals. Each emission i is limited to a rapidity < ]n(M/ku)- Assuming that we can neglect terms ln(fc t i/fc t 2) 
from the rapidity integration, we can write the latter as 

J dy 1 dy 2 Q^y 1 \-\n^Q^y 2 \-\n^Q{R^\y 1 ~y 2 \)=AR\n^-R 2 ^0(R\nC) , (49) 

where £ = kfijktx and we have included the constraint that J ra p(&i, k 2 ) be non-zero. We can then write Eq. ([3]) as 

-dust _ A a2 s ° 2 f 1 d C r ^ f Vi+/+2Cco.* dk tA ( M 2 



-i-^J f J - / \*R* - - R< ] , ,50, 



where we have dropped the O (i?ln£) term of Eq. (1491) . Performing the k t \ integration gives 

„2/-f2 [-X 

C J-* 2tt 



j-clust _ f ^( 



(-2R In — + — ^ lnfl + C 2 + 2( cos <j>) - — ln 2 (l + C 2 + 2( cos < 
Pt 2 J 2 



(51) 



Because J d^ln(l + C 2 + 2C cos <p) = 0, the first term in square brackets vanishes. This was the only term that had a 
NNLL a 2 hxM/pt factor and so at NNLL accuracy J 701 ' 18 * is zero at large R, modulo 1/R corrections associated with 
the 1/R term in Eq. (|47j) . The only element that survives the azimuthal integration in Eq. (|51|) is the second term in 
square brackets, resulting in 

rv 2 C 2 

_ F clust = _ 2 _ i — _R(3. (52) 

This is N 3 LL, so beyond our accuracy. Note, however, that it is enhanced by a factor of R. In the large R limit, 
the separation between partons is limited to be at most 2 In M/pt and thus the R factor is effectively replaced with a 
coefficient of order hxM/pt. Consequently the apparently N 3 LL term of Eq. (|52"j) is "promoted" and becomes a NNLL 
a 2 In M/pt contribution. This is not accounted for in the purely NNLL independent analysis that led to Eq. (|46|) . 
The exact infinite R result can be obtained at order a 2 L by evaluating J" clust with J(fci, k 2 ) = 1, giving 

dust _ 1C « 2 C 2 f 1 d( r d4> dh,i (, M \ , M _ a 2 ^, M , „ M, 



^ciust = 16r s^ / ^ / ^ /v — — rp, ln - lnC ln _^ = _ 4 r| C 3^1n-+0(^ln u -). (53) 

1" JO C J~ir 27T J pt fc t)1 V h,l J «t,l 71" Pt Pt 

Note the agreement of the £3 term here with that derived in Eq. (|4"4")l . It is this contribution that corresponds to the 
C 3 term in Eq. ©. 
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FIG. 3. Comparison of NNLO, NLL+NNLO and NNLL+NNLO results for the jet-veto efficiency for Higgs (left) and Z-boson 
(right) production at 7 TeV. The Higgs plot also includes the result from a POWHEG (revision 1683) [2(J, H3| plus Pythia 
(6.426) [17] 5 EH simulation in which the Higgs-boson pt distribution has been reweighted to match the NNLL+NNLO prediction 
from HqT 2.0 as in [2ll ]. The lower panels show the results normalised to the central NNLL+NNLO curves. 



4. Correlation matrix between 0-jet and inclusive 1-jet cross sections 



As discussed in [21], the prescription that we propose for determining the uncertainties on the 0-jet cross section 
is to treat the uncertainties on the jet-veto efficiency and on the total cross section as uncorrelated. This gives the 
following covariance matrix for the uncertainties of the 0-jet (co-jet) and inclusive 1-jet (<r>i_j t) cross sections: 



eHl + a^l e{l-e)8l 
e{l-e)5l-aH1 (1 - e)H\ 




(54) 



5. Results at 7 TeV 



For completeness, we show in Fig. [3] results for 7 TeV centre of mass energy. The changes relative to the 8 TeV 
results are modest, with very slightly higher efficiencies at 7 TeV. This can be understood because at higher centre 
of mass energy, the PDFs are probed at lower x values, where the scale dependence is steeper, causing the efficiencies 
to drop off more rapidly as one decreases pt,vcto- 



6. R dependence of results 



Figure @] shows the the jet veto efficiency as a function of pt,vcto for several different jet-radius (R) values. Increasing 
the jet radius, more radiation is captured and therefore a jet is more likely to pass the Pt.veto threshold and so be 
vetoed. Consequence the jet- veto efficiency is expected to be lower for larger R values. This is precisely as observed 
in Fig. |U 

Quantitatively, the differences between the R = 0.4 and R — 0.5 results (the values used respectively by ATLAS 
and CMS) are small compared to the uncertainties on the predictions. In contrast, for R = 1 the differences compared 
to the smaller- R results are not negligible. One interesting feature, commented on briefly in the main text, is that for 
the Higgs-boson case, the uncertainties are somewhat smaller for R = 1 than for R = 0.4 and R = 0.5, especially the 
upper part of the uncertainty band. This can be understood with the help of the observation that the upper edge of 
the uncertainty band for the small R values is set by the Q = Mh variant of the resummation (recall that our default 
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FIG. 4. Jet veto efficiency at NNLL+NNLO as a function of pt.veto, comparing several jet-radius values; shown for pp collisions 
at a centre-of-mass energy of 8 TeV, for gluon-fusion Higgs production with Mh = 125 GeV (large mto P limit) and for Z-boson 
production. Uncertainty bands are shown only for 7? = 0.4 and R — 1.0 in order to enhance the clarity of the figure. The 
R = 0.5 uncertainty band is to be found in Fig. [2] The lower panels show the predictions normalised to the central R — 0.5 
results. 



Q is Mh/2). Using Q — Mh increases the size of L. Since the f coml (R) + f cluBt (R) function grows for small R and 
multiplies a 2 s L, a smaller R value magnifies the impact of an increase in Q. 

If, experimentally, one were to consider using larger R values for performing jet vetoes in order to reduce the 
theoretical uncertainties, one concern might be the greater contamination of the jet' s pt from the underlying event 
and pileup. To some extent this could be mitigated by methods such as subtraction |64| . filtering [65| or trimming [66| . 
Note that with subtraction and filtering (when the latter uses two filtering subjets, or more) our jet-veto predictions 
remain unchanged at NNLO and at NNLL accuracy. 



